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We study the magnetic field generation in a neutron star within the model based on the magnetic 
field instability in the nuclear matter owing to the electron-nucleon parity violating interaction. We 
suggest that the growing magnetic field takes the energy from thermal background fermions in the 
neutron star matter. The system of kinetic equations for the spectra of the magnetic helicity density 
and magnetic energy density as well as the chiral imbalance are solved numerically accounting for 
this energy source. We obtain that, for the initial conditions corresponding to a typical neutron 
star, the large scale magnetic field ~ 10^® G is generated during (10"^ — 10®) yr. We suggest that the 
proposed model describes strong magnetic fields observed in magnetars. 
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The most plausible explanation of radiation of soft 
gamma repeaters [ij and anomalous X-ray pulsars 0 is 
the presence of strong magnetic fields B > 10^® G in a 
neutron star (NS). Such highly magnetized NSs are called 
magnetars. Various models, explaining the origin of such 
strong astrophysical magnetic fields, were reviewed in 
Ref. Q- Nevertheless, the issue of the magnetic fields 
generation in magnetars still remains open. 

Recently in Refs. Si we proposed the new model 
for the generation of strong magnetic fields in magne¬ 
tars based on the instability of magnetic fields in dense 
degenerate matter composed of nonrelativistic neutrons 
and ultrarelativistic electrons interacting by parity vio¬ 
lating electroweak forces. The idea that electroweak in¬ 
teraction can induce the magnetic field instability was 
put forward first in Ref. Q. Within our model, bas¬ 
ing on quite natural assumptions about the neutron star 
structure, we could describe the generation of large scale 
magnetic fields, with magnitudes predicted in magnetars, 
during time intervals comparable with magnetars ages. 

Despite the plausibility of the model developed in 
Refs. U it has a significant disadvantage. The in¬ 
stability of a magnetic field, proposed in Refs. ii , is a 
necessary but not a sufficient condition for the magnetic 
field growth. To describe the magnetic field generation 
in magnetars one should indicate the source which feeds 
the magnetic field growth. This issue is addressed in the 
present work. 

In this paper we further develop the model in Refs. [1, 
Q . We start with a brief description of the basic features 
of our model. Then we propose that magnetic fields can 
take the energy from the thermal motion of particles in 
the NS matter. We modify the kinetic equations, derived 
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in Refs. Bi, to account for the magnetic field satura¬ 
tion, and numerically solve them. Finally, we discuss 
our results. In our work we use natural units in which 
h = c = kB = 1- 

Our model is based on the parity violating electroweak 
electron-nucleon interaction (the eN interaction). We 
shall take that the background nuclear matter consists 
of neutrons and protons. This matter is supposed to be 
unpolarized and nonmoving macroscopically. In Ref. Bl 
we derived the averaged effective Lagrangian of the eN 
interaction in the Fermi approximation as 

>Cint = — ^e7° (ElRl + Vr-Pr) V'e: 
VL~K-np(I-4e)](I-20, 

Er = - ^ K - «p(l - 4C)]2^, (1) 

where Gp ~ 1.17 x 10“®GeV“^ is the Fermi constant, 
rin.p are the constant and uniform densities of neutrons 
and protons, t/je is the bispinor electron wave function, 
^ = sin^ ~ 0.23 is the Weinberg parameter, Pl,r = 
( 1 =F 7^)/2 are the chiral projectors, 7 ® = iy^y^y^y^, and 
7^ = ( 7 °, 7) are the Dirac matrices. 

Now let us consider the interaction of ultrarelativistic 
electrons with background matter, described by Eq. o, 
and an external magnetic field B = (0,0, B). The total 
Lagrangian has the form, C = £em + where £em = 
■067^ ipe is the Lagrangian for the interaction 

of an ultrarelativistic electron with the electromagnetic 
field = (0,0,Bx, 0), and e > 0 is the absolute value 
of the electron charge. 

The Dirac equation generated by C was solved in 
Refs. B,0- Using this solution, exactly accounting for 
both the matter interaction and the magnetic field, one 
can compute the induced electric current along the mag- 
neic field Jz = —e('0e7^'0e) + positron contribution, av- 
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eraged using the Fermi-Dirac distribution. This current, 
which is additive to the ohmic current Johm = CTcondE, 
where cTcond is the matter conductivity and E is the elec¬ 
tric field, turns out to be nonzero. If we restore the vec¬ 
tor notations, one gets for the induced electric current 
J = nB. The parameter 11 reads 

n=^(^5 + E5), 

TT 

Ms (mr - ml) , 

Es (Fl - Fr) (2) 


the energy from thermal motion of nucleons and elec¬ 
trons, which NS is composed of. For this purpose we 
shall calculate the temperature corrections to the energy 
density of degenerate fermions in NS as a possible source 
for the growth of the magnetic field. 

We shall start with the electron component of NS mat¬ 
ter. Using the expansion of the integral in Ref. Q, 


fds_ M _ 

X exp[(s - M)/r] + 1 

= l m<ie + -T^ 


d/(£) 


+ 0{T^), 

e—(i 


(5) 


where Oem = e^/47r Ri 7.3 x 10“^ is the fine structure 
constant. Note that, since we consider ultrarelativistic 
electrons, we can assume that right and left chiral projec¬ 
tions of the electron-positron field behave independently 
and possess different chemical potentials: /iR and /in- To 
obtain Eq. m we assume that ^ Up inside NS. 

Using Eq. m, in Ref. Q we derived the system of 
kinetic equations for the spectra of the magnetic helicity 
densiy h{k,t) and magnetic energy density p^(k,t) as 
well as the chiral imbalance p^{t) in the form. 


dh{k, t) 
dt 

dpB{k,t) 

dt 

dM5(0 

dt 


2fc2 


^cond 

2fc2 


<^cond 


h{k,t) + ( ) PBik,t), 

\ <^cond / 

PB{k,t) + ( ^ ) k'^h{k,t), 

\ ^cond / 

J dk [k'^h{k,t) — 2IlpB{k,t)^ 


- T/Ms, 


(3) 


where p is the chemical potential of electrons in NS, F/ = 
daemm-e/STTCTcond is the chuality flip rate in the electron- 
proton (ep) collisions, and me is the electron mass. Note 
that the chirality flipping term in Eq. m should contain 
/is since the equilibrium in the system of right and left 
electrons is achieved when /tr = pb- 

The total magnetic helicity H and the magnetic field 
strength B can be found on the basis of h{k,t) and 
PB{k,t) as 


H{t) = j d^x{A • B) = F j h{k,t)dk, 

dkpBik,t), 


(4) 


where F is the normalization volume and the integration 
is over all the range of the wave number k variation. It 
should be mentioned that in Eqs. m and O we assume 
the isotropic spectra. 

In Ref. 0 we found that the model described by Eq. m 
reveals the potential growth of the seed magnetic field 
Bq = 10^^ G up to R > 10^^ G, i.e. the strengths pre¬ 
dicted in magnetars. However, the energy source feeding 
the magnetic field growth was not specified in Ref. 

We demonstrate below that the magnetic field can take 


one gets for the energy density 

f _ P __ , c 

J (27r)3 exp[(p -/i)/r]-b I 

P-o Spe = (6) 

and the number density 



d^p _1_ 

(27r)3 exp[{p- p)/T] + 1 
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of degenerate ultrarelativistic electrons including tem¬ 
perature corrections. In Eqs. (jb)) and 0 we keep only 
the leading terms in the temperature T. To derive 
Eqs. ([6]) and d?]) we neglect the magnetic fields correc¬ 
tion to pe and He, studied in Ref. since eB <C p^ for 
B = (10^^ — 10^^) G we consider here. 

One can see in Eqs. ([6]) and 0 that the mean energy of 
a thermal electron (ee)r = SpejSne = 3p/2 exceeds the 
Fermi level p. The cooling of such electrons proceeds in¬ 
dependently of the main contribution in degenerate elec¬ 
tron gas with 0 < Ee < p since both the energy density 
of electrons and their number density are proportional to 
T^. This cooling does not violate the Pauli principle for 
them either. That is why the decreasing of the tempera¬ 
ture of such thermal electrons can feed the magnetic field 
growth. 

Now let us consider degenerate nonrelativistic nucleons 
TV, i.e. neutrons N = n and protons N = p, as the energy 
source for the magnetic field growth. These particles have 
the Fermi energy p^ = p^p^/2Mfq T, where pfn is the 
nucleons Fermi momentum and Mjq is the nucleon mass. 
Analogously to Eqs. ([6]) and 0, as well as using Eq. 0, 
we get the energy and number densities for degenerate 
nucleons, including thermal corrections, as 


Pn =2 


d^p 


Pno = 


P% 


(27r)3 exp[(e - pn)/T] -b I 
T'^MmPFn 


= Pno + SpN, 




, SpN = 


( 8 ) 
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and 


riN 

I^NO 


=2 


d^p 


1 


(27r)3 exp[(£: - pn)/T] + 1 
Stt^ ’ ^ 6pfn 


P% 


nNo + dnjv, 


(9) 


where we account for the energy-momentum relation for 
a nonrelarivistic nucleon e = p^/2Mjs[ and again keep 
only the leading terms in T. 

Basing on Eqs. and one obtains the mean 
energy of thermal nucleons (eAr)T = Spn/Sun = 
3pp^/2Mn, which is above the Fermi surface: {ejsf)T > 
Pn- Hence, like electrons, these nucleons can transfer 
their thermal energy to the magnetic field in their cool¬ 
ing without violation of the Pauli principle. 

Summing up the thermal energy density corrections 
of electrons, protons, and nucleons, we can define the 
equipartition magnetic field strength i?eq as 

=Spe + Spp + Spn 

■ MnPFr, + MpPF^ 

2 

Accounting for Mn r; 940 MeV, ~ 

339 MeV for the NS density n„ = 0.18 fm“^, and pFp ~ 
p = 125 MeV for the electron density rie = 9 x 10^® cm“^, 
one gets that the neutron contribution to is the great¬ 
est one. We can consider the quantity px = -Sgq/2 in 
Eq. (HOI) as the inexhaustible energy source requiring that 
Hgq B^. Thus we do not violate the total energy 
conservation for the extended system which includes the 
background matter and the magnetic field. Since the val¬ 
ues of rin.e are typical for NS we shall later use them in 
the numerical simulations. 

In Refs. S0 we simulated magnetic fields in magne- 
tars solving Eq. (jS]) and using Bq = 10^^ G as the initial 
magnetic field. Assuming T = 10® K and the above pa¬ 
rameters of the NS matter, we get that Bq <C Hgq- How¬ 
ever, if H = 10^^ G, one obtains that B^ ^ B^^. Thus 
strong magnetic fields, predicted in Refs. ii , will in¬ 
fluence the background matter in NS. 

To avoid a back reaction on matter from such a strong 
magnetic field we should modify Eq. As known from 
the solar dynamo theory Q, one can avoid the infinite 
growth of the magnetic field by quenching of the dynamo 
a-parameter. Thus we can introduce the quenching of 
the parameter H, given in Eq. @ , as 



H^ 


H 


( 11 ) 


where and Rgq can be found in Eqs. O and m- 
Again referring to the solar dynamo theory, Heq is equiv¬ 
alent to Bq ^ 1 kG, which is the magnetic field strength 
in a solar spot. Now the excessive growth of the mag¬ 
netic field is eliminated from our model since it is the 


parameter H which is responsible for the magnetic field 
instability. 

To analyze the magnetic field generation in a magnetar 
on the basis of Eq. ([3]) we should adopt an appropriate 
initial condition. The detailed discussion of the initial 
condition is provided in Ref. Here we just make a 
few comments on it. 

We shall consider the evolution of a thermally relaxed 
NS at t > to, where to = 10^ yr. As obtained in Ref. 0, 
at ^0 < ^ ^ 10® yr, NS cools down by the neutrino emis¬ 
sion in modified Urea processes. The time dependence of 
the NS temperature obeys the differential equation [HI, 

dT(t) r(t) 

dt (nT-2)t’ ^ ’ 

where the index np = 8 for modified Urea processes. Us¬ 
ing Eq. (HU and the results of Ref. Q, one gets that 
the NS temperature and the NS conductivity will de¬ 
pend on time as = T^F and CTgond = o'o/F, where 
F = (t/to)“i/®, To = 10® K, and cto = 2.7 x 10® MeV 
is given by the electron (or proton) density Ue = rip = 
9 X 10®® cm"®. 

We shall study the generation of the magnetic field 
without specifying its direction, which can be random. 
Moreover we suggest that a seed magnetic field appears 
due to a turbulence which can be of a hydrodynamic ori¬ 
gin. In this case one can choose the initial Kolmogorov 
spectrum for the magnetic energy density PB(fc,^o) = 
Cfc“®/® fo . Here we correct the initial spectrum chosen 
in Ref. The constant C can be found from Eq. o 
setting B(to) = Bg = 10^^ G, which is a seed field typical 
for a young pulsar. The wave number runs in the inter¬ 
val /Cniin < A: < A:niax, where /Cmin — -^NS ~ ^ ^ 10 ^V, 

i?NS = 10 km is the NS radius, and Ab is 

the free parameter specifying the scale of the magnetic 
field generated. 

The initial spectrum of the helicity density can be cho¬ 
sen as h{k,tg) = 2qp^{k,tg)/k, where the parameter 
0 < g < 1 defines the initial helicity. The case g = 0 
corresponds to the initially non-helical field and g = 1 
to the magnetic field with a maximal helicity. Therefore, 
besides magnetic fields we can also study the generation 
of the magnetic helicity in our model. 

The initial value of the chiral imbalance can be taken 
as pg{tg) = IMeV. Note that pg{tg) ^ 0 is generated 
in direct Urea processes at the early stages of the NS 
evolution at f < to- The energy scale of these processes 
is governed by the mass difference between a neutron and 
a proton: — Mp ^ 1 MeV. This fact substantiates our 

choice of pgitg). 

At the first glance one can imagine that, for the chosen 
parameters, the contribution of the electroweak interac¬ 
tions ^ Vs to Eq. ([3]) is negligible compared to the elec¬ 
trodynamic contribution ^ pg. As shown in Refs. S0, 
almost any initial pz(tg) 0 tends to zero very rapidly 
because of the high rate of ep collisions, whereas 14 is a 
steady source for the growth of the magnetic helicity and 
the magnetic energy density. Moreover, the electroweak 
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term is not affected by ep collisions since V 5 depends 
on the difference of the interaction potentials of left and 
right electrons with background matter, which are con¬ 
stant parameters of the model [see Eqs. (HD and @] , un¬ 
like /Ts, which is a dynamic variable. 

We also mention the recent Ref. [l^, where another 
steady source for the magnetic field instability, different 
from V 5 , was used to explain strong magnetic fields in 
magnetars. It is based on the generation of the chiral 
imbalance in direct and modified Urea processes, + 
p ^ n + I'eL and + p + N —>■ n-|- i/eh + N, which 
are not in the equilibrium with inverse reactions. This 
situation can happen during ^ 10 s after the onset of the 
supernova collapse outside the neutrinosphere. 

Short scale, Ab < 1 cm, magnetic fields with the 
strength B < 10^'* G were demonstrated in Ref. [T^ to 
be generated in this situation. However, as shown in 
Ref. [ 1 ^, short scale chaotic magnetic fields in a super¬ 
nova explosion are subject to the reconnection with the 
typical time of several seconds. This time scale is compa¬ 
rable with the time interval for magnetic field generation 
in Ref. [3. Thus, magnetic fields predicted in Ref. [l^ 
will transform effectively into heat because of the mag¬ 
netic reconnection. 


Below we present the results of numerical solution of 
Eq. ([3]) accounting for Eq. (fTTl) and the chosen initial 
conditions. In Eig. [T]one can see the growth of magnetic 
fields of different length scales and initial helicities. We 
study the two main minimal scales: = I km in 

Figs. 1(a) and |l(b)| as well as a|^™^ = 100 m in Figs. 1(c) 


and 1(d) Thus we predict the generation of strong large 
scale magnetic fields. 

To compare the behavior of magnetic fields in the 
present work with that in Ref. Q, in Figs. 1(a) and |l(c)| 
we also show the results of the numerical solution of 
Eq. ([3]) without quenching in Eq. (fTTl) . One can see 
that unquenched magnetic fields, shown by blue l ines, 
slow down the growth rate after ^ 10® yr in Fig. 1(a) 
and ^ 10^ yr in Fig. 1(c) but continue growing |3J. 
On the contrary, the quenched magnetic fields, shown 
by red lines, are saturated. For both A^'"^ we start 
with Bq = 10^^ G and magnetic fields reach the satu¬ 
rated value Bsat ^ 10^® G. For example, in Fig. |l(b)[ 
Bsat ~ 1-1 X 10®® G. This Rgat is close to magnetic fields 
observed in magnetars (T^ . 

Magnetic field s in F ig. 1(a) grow up to Bsat for t > 
10® yr and in Fig. 1(c) for t > 1C)‘® yr. These time intervals 
are comparable with the ages of young magnetars [13- 
Note that the smaller the scale of the magnetic field is, 
the faster this magnetic field grows and the stronger Bsat 
is. One gets from Eq. (|3]) that pB ^ where A is 

the typical vector potential. Hence, a bigger fcmax corre¬ 
sponds to stronger Hsat- 

We also analyze the evolution of magnetic fields with 
different initial helicities. One can see in Figs. |l(b)| 
and |l(d)| that there is a difference in the behavior of 
magnetic fields for initially non-helical (solid lines) and 


maximally helical (dashed lines) fields for relatively small 
evolution times. At later times this difference is washed 
out; cf. Figs. |l(a)| and 1(c) It means that, in frames of 
our model, we can generate both strong magnetic fields 
and the magnetic helicity. 

Note that the behavior of quenched and unquenched 
magnetic fields is almost indistinguishable at small evolu¬ 
tion times. Indeed, if t <C tsat, where fgat = (10^—10®) yr 
is the saturation time depending on the scale of the mag¬ 
netic field, then B <C i?eq in Eq. dinD- Thus in this 
time interval it is sufficient to consider the evolution of 
quenched magnetic fields, which is shown in Figs. |l(b)| 
and 1 1(d)] 

Gomparing the results of Ref. [3 w ith the evolution of 
magnetic fields in Figs |l(b)] and |l(d) one can notice that 
in the present work magnetic fields grow several times 
slower. This discrepancy can be attributed to the fact 
that now we use the correct initial Kolmogorov’s spec¬ 
trum PB{k,to) = Ck'^'^ with vb = —5/3 vs. vb = 1/3 in 
Ref. [ 3 . Indeed, since C ^ + 1, the greater vb is, the 

faster pB{k,t) will grow. 

Along with growing magnetic fields, shown in Fig. [l] 
it is important to consider the evolution of the magnetic 
helicity density h{t) = H{t)/V to illustrate its generation 
in a magnetar. In Fig. [3] we demonstrate how the mag¬ 
netic helicity grows in our model. We consider the cases 
of initially helical and nonhelical magnetic fields as well 
as quenched and unquenched parameter H in Eq. m to 
compare our results with those in Ref. [3 . One can see in 
Figs. |2(b')] and [2(d)] that the difference in the evolution of 
initially helical and nonhelical magnetic fields is impor¬ 
tant only at early evolution times. Later this difference 
is washed out; cf. Figs. 2(a) and |2(b)[ Therefore we ex¬ 
tend the result of Ref. [3j, that the magnetic helicity can 
be generated in our model, to the case of the quenched 
parameter H in Eq. (ED- 


In conclusion we mention that we have further devel¬ 
oped the model, recently proposed in Refs. Si , for 
the magnetic fields generation in magnetars. We have 
improved our approach pointing out that the magnetic 
field, growing owing to the instability caused by the par¬ 
ity violating eN interaction, can take the energy mostly 
from thermal neutrons, as well as electrons and protons, 
which are present in the NS matter. 


We have started with the evaluation of thermal cor¬ 
rections to the number densities and the energy densities 
of background fermions in NS. We have shown that, by 
cooling, these particles can pass their thermal energy to 
the magnetic field without violating the Pauli principle. 
Then, in the analogy with the solar dynamo, we have 
generalized the kinetic equations, derived in Ref. i,by 
quenching of the parameter H; cf. Eq. ED- This proce¬ 
dure allowed us to treat background fermions as the large 
energy reservoir feeding the magnetic field. Moreover we 
have avoided the infinite growth of the magnetic field. 

We have numerically solved the system of kinetic 
Eqs. with the modified parameter H. For the ini¬ 
tial conditions corresponding to a typical NS {nn,e and 
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FIG. 1. (color online). The growth of quenched and unquenched magnetic fields in magnetars versus t — to for different fcmax, 
corresponding to different length scales Ab . Red and blue lines in panels (a) and (c) are the solutions of Eq. (|3} for quenched 
and unquenched 11 in Eq. m respectively. Solid lines correspond to initially nonhelical fields with q = 0 and dashed lines 
to maximally helical fields with q = 1. (a) The magnetic field evolution for fcmax = 2 x 10“^° eV or = 1km. (b) The 

behavior of the quenched magnetic field with the parameters as in panel (a) for shorter evolution time to < t < 5 x 10® yr. 
(c) The magnetic field growth for fcmax = 2 x 10“® eV or A^‘"^ = 100 m. (d) The behavior of the quenched magnetic field with 
the parameters as in panel (c) for shorter evolution time to < t < 5 x 10® yr. 


Bq), we have obtained the growth of the seed magnetic 
field by three orders of magnitude to Bgat ~ 10®^ G. Al¬ 
though this value of Bsat is smaller than that obtained 
in Refs. Si, this Rs at. i s close to the magnetic field 
predicted in magnetars |lq . 

The time of the magnetic field growth to Bsat is Gat = 
(10"^ — 10®) yr depending on the scale of the magnetic 
field. We have analyzed the two scales of the magnetic 
field in the range Ab = (10® — 10®) m, i.e. we predict the 


generation of large scale magnetic fields. Comparing the 
obtained results for tgat with the ages of magnetars [ij . 
one concludes that our model is a quite plausible expla¬ 
nation of magnetic fields in magnetars. 
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versity and to RFBR (research project No. 15-02-00293) 
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FIG. 2. (color online). The evolution of the magnetic helicity in magnetars versus t — to for quenched and unquenched parameter 
n as well as for different fcmax, corresponding to different length scales Ab. Red and blue lines in panels (a) and (c) are the 
solutions of Eq. ® for quenched and unquenched 11 in Eq. HU respectively. Solid lines correspond to initially nonhelical 
fields with q = 0 and dashed lines to maximally helical helds with q — 1. (a) The magnetic helicity density evolution for 
fcmax = 2 X 10“^° eV or = 1 km. (b) The behavior of the magnetic helicity density for the quenched 11 with the parameters 

as in panel (a) for shorter evolution time to < t < 5 x 10® yr. (c) The magnetic helicity density growth for fcmax = 2 x 10~® eV 
or Ag™‘"^ = 100 m. (d) The behavior of the magnetic helicity density for the quenched 11 with the parameters as in panel (c) 
for shorter evolution time to < t < 5 x 10^ yr. 
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